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ABSTRACT 

We present binary galaxy merger simulations with varying mass ratios and different progenitor 
morphologies. The simulations include mergers of gas-rich disks (Sp-Sp), of early- type galaxies and 
disks (E-Sp, mixed mergers), and mergers of early-type galaxies (E-E, dry mergers). We follow the 
dynamics of gas, stars and dark matter, and include radiative cooling, star formation and black hole 
(BH) accretion and the associated feedback processes as in Springel et al. We study the mass assembly 
of the BHs and discuss technical issues of the implemented model in detail. For Sp-Sp mergers, the 
peak star formation rate and BH accretion rate decrease and the growth timescales of the central black 
holes and newly formed stars increase with higher progenitor mass ratios. The peak BH accretion 
rate typically occurs shortly after the time of BH merging for low progenitor mass ratios (e.g. 3:1 
and lower), whereas for higher progenitor mass ratios there is no clear correlation between the peak 
BH accretion rate and BH merging time. The termination of star formation by BH feedback in disk 
mergers is significantly less important for higher progenitor mass ratios (e.g. 3:1 and higher). In 
addition, the inclusion of BH feedback suppresses efficiently star formation in dry E-E mergers and 
mixed E-Sp mergers. All merger remnants, independent of their progenitors, follow the observed 
relations between the central BH mass and the stellar velocity dispersion {Mbh — c), the bulge mass 
(Mbh — M^,) and the bulge binding energy {Mbk — M^,a'^), with the dominant source of scatter arising 
from variations in the initial gas mass fraction. The normalizations for all relations and the simulated 
slope of the Mbh — c and Mbh — ^ relations are in good agreement with the observations, whereas 
the simulated slope of the Mbh — Af* relation is slightly steeper compared to the observations. This 
indicates that the simple BH feedback model self-regulates the mass growth of the BHs in accordance 
with the observed relations for a very wide range of merger progenitors. 

Subject headings: galaxies: interaction- galaxies: active galaxies: evolution - galaxies: formation — 
methods: numerical 



1. INTRODUCTION 

Observations in the last decade have shown that 
supermassive black holes (BH) reside in the cen- 
ters of most if not al l gala xi es with spheroids (e.g . 
iKormendv &: Richstond 119951 : iRichstone eFall Il998h . 
The properties of these black holes and their host galaxies 
are correlated in a number of ways, e.g. the correlation 
of the BH mass and the bulge veloc ity dispersion as ex- 
press ed in the Mbh — o" r e lation f e.g. iFerrarese fc Merrittj 
120001: iGebhardt et"al] 120001 : iTremaine et all I2002D . 
the relation between BH mass and the bulge stel- 
lar mass, Mbh — Mb ,, ]^^ fe.g. IMagorrian et al.l 119981 : 
IMarconi fc HuntI 120031: IHaring fc Rixll2004D and the re- 
lation between BH mass and the bulge binding en- 
ergy, Mbh - Mhnigscr^ fe .g. lAller fc Richstond 120071 : 
iBarwav fc Kem bhavi'2007^. These observed relations in- 
dicate that supermassive BHs and the stellar spheroids 
of their host galaxies do not form independently of each 
other. 

There exists several theoretical explanations for the 
observed correlation between galaxy properties and the 
black hole mass, including physical mec hanisms such as 
visco us disk accretion and star forma tion (iBurkert fc SiT3 



and their asso ciated interm ediate-mass BHs to the 
galac tic centers (jBournaud e t al. 2007a; El megreen et al.l 
|2008[) . The majority of the models rely on some 
form of self-regulated BH mass growth, in which gas 
is fed to the central black hole until the black hole 
releases sufficient energy to unbind the gas and blow 
it away in momentum - or pressure-driven winds fe.g . 
Ciotti fc O strikcr 1997: Silk fc Rees '199^: iFabianI [l999l : 



Murray et a l. 2005; Ciotti fc Ost riker 200 



[200l . gas or dark matt er collapse (I Adams et al.ll2001f ). 
dissipationless merging (jCiotti fc van Albadall200l[r and 

the merging and migration of massive stellar clumps 
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The first theoretical studies of the co-evolution of black 
holes and galaxies were carried o ut using sem i-analytic 
modeling fSAM) techni q ues fe .g . fCattaneojj t al. 199!|; 
Kauffmann fc Haehnclt 20 001: IWyi thc fc Locb 200|: 
VoloQtcri et al. 20 03; Grau ato et al.l 12004: Croton et al] 
20061 iBower et all [20061) . In SAMs the evolution of the 
dark matter halos is followed with merger trees and the 
evolution of the galaxies and their embedded black holes 
is modeled with simple physical parameterizations. Most 
SAMs are based on the assumption that quasar activity is 
triggered by major galaxy mergers of gas-rich galaxies, in 
which large quantities of gas flow to the cente rs of galax- 
les (e .g. fe ar nes & Hernquist 1996; Kazant zidis et all 
120051 iNaab et al.ii2006at : iHopkins et al.i i2008b') produc- 
ing starbursts a nd leading to rapid black hole growth 
(|HernauistlU"989l ). During the mergers the dominant dy- 
nami c al process is the violent relaxation (jLvnden-Belll 
Il967t lArad fc JohanssonI l2005f ) of the stellar and dark 
matter components. 

This picture is supported by observations of local merg- 
ers that show large central gas concentrations, accompa- 
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nied by strong starburst and AGN activity (|Genzel et al.l 
Il998l ). indicating BH growth. These so called ul- 
traluminous infrared galaxies (ULIRGs) mostly origi- 
nate from major mergers, typically wi th mass ratios 
below 3:1 (e.g. Dasv ra et al.l I20d6al l3: IVaisanen et al.l 
[2008,). The ULI RGs could eventually evolve into 
spheroidal galaxies (iHibbard fc YunI 119991: iGenzel et al.l 
120011: iRothberg fc JosephI 120061: iDasvra et al .' '2006bl). 
In addition, nearby quasa rs fe.g. iMcLure e t al. 1999; 
ICanalizo fc StocktonI l2001h are preferentially found in 
tidally disturbed objects, thus corroborating the impor- 
tant role of galaxy interactions and mergers for major 
black hole growth. Furthermore, an increasing number 
of recent infrared and X-ray surveys have completed the 
census of BH activity and mass growth by revealing a 
significant populatio n of highly obscured AGNs at red- 
shifts above z \ (I Johansson et al.l 120041 : iHouck et al.l 
[20051 iBareer et al.ll2005D . 

A first, simple self-consistent method for includ- 
ing black hole growth and thermal feedback coupling 
to the surrounding gas in sim ulations of binary dis k 
galaxy mergers was developed bv lSpringel et al.l (|2005bD . 
Based on this description hydrodynamical simulations 
of BH growth in galaxy mergers w ere used to re- 
produ ce the local Mbh — tr relation (|Di Matteo et al.l 
120051). study the reds hift evolution of this relation 
([Robertson et al.ll2006cf ). provide an explanation of th e 
red colors of massive ellipticals ([Springel et al.l [2C)05aD . 
and describe the fundam ental plane of elliptical galaxies 
([Robe rtson ct al." 2006b). Using this black hole merger 
sample Hopkins ct al. (2005, 2006) proposed a compre- 
hensive picture of a unified merger-driven model for the 
origin of starbursts, quasars and their relation to galaxy 
spheroid formati on. Recently, simulati ons of both indi- 
vidual galaxies (fKhalatyan et al.' '2008*) and larg er cos- 
molo gical volumes (Sijacki ct al. 2007; Di Mattco ~et al.l 
l2008f l with BH feedback starting directly from initial con- 
ditions appropriate for the ACDM cosmology have been 
performed. Both cosmological studies assumed that seed 
black holes were present at high redshifts in forming halos 
and then followed their evoluti on self-cons istently with 
hydrodynami cal simulations . [Pi Matteo et al.l (|2008D 
employed the lSpringel et al.l (| 2005bl) black hole feedback 
prescription, whereas ISiiacki et al.l ([2007f ) used in ad- 
dition a predominantly mechanical feedback model for 
black holes in galaxy clusters with low accretion rates, in 
which AGN-d riven bubbles are i njecte d into a gaseous 
environment (jSi iacki fc Sp ringell l2006f ) . The study by 
iDi Matteo et al . (2005) showed that the Mbh - o rela- 
tion can be reproduced for equa l-mass disk galaxy m erg- 
ers. The cosmological s tudies of lSiiacki et al.l ([2007[ ) and 
iDi Matteo et al.l ([2008D naturally include mergers with a 
range in mass ratios and morphologies, and both studies 
are able to find adequate fits to the observed Mbh — cr 
and Mbh — M* relations at both high and low redshifts, 
where M, is defined as the total stellar mass within the 
effective radius. However, even the highest resolution 
cosmological simulations resolve spatial scales that arc 
an order of magnitude larger than the typical spatial 
scales resolvable in high-resolution merger simulations 
(e ~ 0.1/i~^kpc) and it is not obvious that low-resolution 
simulations that only resolve scales corresponding to the 
effective radii for typical galaxies are able to capture the 
essential physics related to BH feedback. 



In a cosmological context unequal-mass mergers are 
more frequent than equal-mass mer gers and the progeni- 
tors will have varyin g morphologies ([Khochfar fc BurkertI 
120011: iKhochfar fc"S ilk 2006). It has been shown that 
the mass ratio of disk mergers has a significant impact 
on the remnant properties. Equal-mass disk mergers 
predominantly result in slowly rotating systems with 
boxy or disky isophotes, whereas unequal-mass merg- 
ers lead to rotating, i ntrinsically an i sotropic spheroids 



with disky 
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is ophotes ([B ar 
n[200lt [B arnes 



arne^ '19981 iNaab et al 



2 002: Naab fc BurkertI 



199£ 



200S 



Bournaud et al.l 1200 5: Burk ert fc NaabI 120051 : ICox et all 



20061 : iBurkert et an[2007i) . The mass ratio and the actual 



gas content have a significant influence on the final stel- 
lar orbit st ructure and in ternal kinemati cs of the merger 
remnants (iBarnes fc Hern auist 1996.; Jesseit et al.ll2005l 
120071: iThomas et al.ll2007D . Apparently, binary mergers 
of disk galaxies with varying mass ratios are one rea- 
sonable mechanism for the formation of intermediate- 
and low-mass ellipticals, whereas more massive slowly 
rotating ellipticals are probably formed in multiple merg- 
ers a t high redshifts ([Naab et all 120071: IBurkert etall 
|2007| ). A potential assembly mechanism for massive 
boxy and slowly-rotating elliptical galaxies at low red- 
shift is mergi ng of gas-poor early-type galaxies (dry 
mergers) (Na ab et al.l l2006b( l. There is now mount- 
ing observational e vidence for dry galaxy rnergers both 



in ga laxy clusters ([van Dokkum et alii 199 91: 'Tra n et al 
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2005D and in the field ([van DokkumI 12005: Bell et al 



20061 ). Simulations of gas-free mergers have been shown 



to preserv e the fun damental plane of elliptical galaxies 
([Nipoti et a l. 2003: iBovlan-Kolchin et all l2005f ) . Fur- 



thermore [Robertson et al. I ([2006bf ) showed that E-E 
mergers including BH feedback physics do not strongly 
alter the original Mbh — f relation, although they might 
be the source for some of the observed scatter in the re- 
lation. Mixed merging, in which a gas-poor early-type 
galaxy merges with a gas-rich spiral disk galaxy might 
also form elliptical galaxies with distinct properties . 
They were shown by Khochfar fc Burkert[ ([2001 [2005h : 
Naab ct al. (2006b); Hopkins ct al. (20083) to dominate 
the merger budget at intermediate masses between the 
low-mass ellipticals Sp-Sp mergers and the highest mass 
E-E mergers. Thus it is of utmost importance to study 
whether the current BH feedback model that reproduces 
the Mbh — cr and A/bh — M, relations for equal-mass 
disk mergers, is also able to reproduce the observed re- 
lations in galaxy mergers with varying mass ratios and 
progenitor morphologies. 

In this paper we present simulations of unequal-mass 
disk mergers, mixed E-Sp mergers and E-E mergers, in- 
cluding BH feedback. Our aim is to study the observed 
local A/bh — CT and A/bh — Af, relations and see whether 
they hold for our merger sample with varying mass ra- 
tios and different progenitor morphologies. We also study 
in detail the effect of varying the merger mass ratio on 
the peak BH accretion rate and final BH mass. In ad- 
dition, we analyse the termination of star formation by 
BH feedback in unequal-mass disk mergers, mixed E-Sp 
mergers and dry E-E mergers comp ared to equal-mass 
disk mergers ([Springel et al.l [2005a[ ). Furthermore, we 
discuss in deta il the modeling of black hole mergers in the 
ISpringel et al.l (i2005b) feedback prescription. The merg- 
ing of black holes is not resolved in the numerical simu- 
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lations and therefore an effective model for BH merging 
needs to be adopted. We show that the assumption of 
rapid merging of the black holes in equal- and unequal- 
mass galaxy mergers is critical for reproducing the final 
BH mass that fulfil ls the observed relations using the self- 
regulated .SjjringeLiE^n |200 5b) BH feedback model. 

This paper is structured as follows. In ^we describe 
the simulation code and discuss the BH feedback and 
merger prescriptions. The galaxy model and merger se- 
tups, together with the adopted parameter settings, are 
discussed in [J3l Here we also simulate a set of isolated 
disk galaxies and perform a detailed analysis of the effect 
of the merger prescription and resolution on the growth 
and final mass of the BHs. In [2] we present the re- 
sults and implications of our Sp-Sp, mixed E-Sp and dry 
E-E simulation sets to the origin of the A/bh — cr, the 
Mbh — Af* and the A/bh — Af^cr^ relations. In addition 
we perform here a detailed study of BH and star forma- 
tion activity as a function of merger mass ratio. Finally, 
we summarize and discuss our findings in SJS] 

2. METHODOLOGY 
2.1. Numerical code 

We perform the simul ations using the parallel 
TreeSPH-code GADGET-2 (|Sr)ringe]|[2005h . The code 
follows the gas dynamics using the Lagrangian Smoothed 
Particle Hydrodynamics (SPH) (e.g. iMo naghan 1992) 
technique formulated in such a way that energy and 
entro py is manif estly conser v ed (|Springel fc Hernquis^ 
I2002D . Following iKatz et al.l (|1996f) the radiative cool- 
ing for a primordial mixture of hydrogen and helium to- 
gether with a sp atially uniform time-in dependent local 
UV background (|Haardt &: Madaul[l996l ) is included. 

We include star formation and the associated super- 
nova feedback, but exclude supernova-driven galactic 
winds, following the sub-resolution mult iphase model de- 
veloped by Springel & Hcrnquist| ()2003[ ). In this model 
the ISM is treated as a two-phase medium in which cold 
clouds are embedded in a tenous hot gas at pressure equi- 
librium. Stars form from the cold clouds in regions were 
p > Pth with the shortlived stars supplying an energy 
of 10^"'^ ergs to the surrounding gas by supernovae. The 
threshold density, pth, is determined self-consistently in 
the model by requiring that the equation-of-state (EOS) 
is continuous at the onset of star formation. 

We have implemen ted all the basic aspects of the 
ISpringel et al.l ()2005bD effective BH feedback model into 
our version of GADGET-2. We summarize the main as- 
pects of the BH feedback model in Appendix |Al In Ap- 
pendixlBlwe perform a direct comparison of our feedback 
implementation with the original model and find that the 
two models are in very good agreement. 

2.2. Modeling black hole mergers 

In the final stages of a galaxy merger the dark mat- 
ter halos of the galaxies merge and in the center a sin- 
gle stellar system is formed from the coalescence of the 
two stellar components. During this process a substan- 
tial fraction of the gaseous component is funneled to the 
center where it can drive strong starbursts and AGN ac- 
tivity. Presumably during this process the two central 
supermassive BHs form a binary, that hardens and even- 
tually leads to a coalescence of the BHs. However, it 



is not clear how long this process would last and what 
the dominant mechanism for the hardening of the binary 
would be. Possible mechanisms includ e hardening of the 
BH b i nary through stellar-dynamica l (jMakino fc Funatol 
12004 iMerritt fc Milosavlievic'2005') and hydrodynami- 
cal (jEscala et al. 2004, ) processes. At the final stages 
of the BH merger {d < O.Olpc) gravitational radiation 
will dominate angular momentum a nd energy losses and 
cause the two BHs to coalesce (e.g. iFlanagan fc Hughe^ 
Il998f) . However for unequal-mass BH pairs, gravitational 
wave emission also removes net linear momentum and 
can impart a recoil velocity to the center of mass of the 
system. This recoil velocity is typically of the order of 
'^CM ~ 100 — 400 kms^^ and may play a role in removing 
relatively low-mass supermassive BHs from d warf galax- 
ies at moderately to high redshifts (Madau fc Quataertl 
[200l[MCTritt et al.l[200l . 

In the simulations presented in this paper we lack the 
required resolution to follow the hardening and eventual 
merging of the BHs. Using current numerical techniques 
it is not possible to accurately resolve the dynamical 
evolution of a galactic BH pair on the relevant scales, 
although alterna tive techniques su ch as the particle- 
splitting method (iMaver et al.ll2007f) and a daptive mesh 
refinement simulations ( Levine et al.l 120081 ) might allow 
one to follow the BH dynamics to somewhat smaller sep- 
arations compared to traditional methods. Instead, as 
was the case for the BH accretion, we use an effective 
model that ensures rapid merging of the BHs once they 
come sufficiently close to each other. There is circum- 
stantial evidence suggesting that the merging of BHs is 
efficient. If this was not the case, it should be common 
for a third black hole to be brought in by merger, in 
which case the BH binary and thir d BH could lead t o 
sling-shot ejection of all three BHs (jSaslaw et al.lll97^ . 
This process would seriously hamper the growth of su- 
permassive black holes. Thus the very existence of super- 
massive black holes together with the rel atively low num- 
ber of known supermassive BH bina ries ([Sillanpaa et al.l 
[19881: iMCTritt fc MilosavIie^ [2005 11 can be taken as an 
indication that BH binaries typically merge efficiently. 
Upcoming gravitational wave experiments, such as LISA 
will bring more clarity to the duratio n and true number 
of BH mergers in the local Universe (jHugh cs 2002)^ 

In this study we adopt the criterion of i Springel et al.l 
(2005b) for BH merging, which allows for rapid BH merg- 
ing. Two BHs are assumed to merge instantly if they 
come within the smoothing length of each other and if 
their relative velocity at this time is below the local sound 
speed, as determined from the nearest neighbors within 
the SPH kernel. The momentum is conserved when gas is 
absorbed by the BH sink particle, thus an accreting BH 
moving in a gas-rich environment will experience a fric- 
tion force that reduces the relative motion of the BH with 
respect to the surrounding gas. However, we found that 
for unequal-mass mergers this description was not always 
adequate in ensuring rapid merging of the black holes. 
Typically the black holes had too large relative velocities 
to allow for merging and thus an alternative method had 
to be employed to ensure the merging of the BHs. After 
discussions with Volker Springel (Springel, private com- 
munication) we decided to adopt a model, in which the 
BH is repositioned at every time step to the minimum of 
the local potential. This minimum is evaluated by find- 
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ing the gas particle within the SPH smoothing length of 
the BH which has the minimum gravitational potential. 
The BH is then repositioned to the location of this min- 
imum at every timestep. This method effectively glues 
the BH to the center of the galaxy. Using this description 
the BH will remain in the merging galaxies as long as the 
galaxy core structure is intact. After disruption the two 
BHs will rapidly slide to the center of the parent galaxy 
and merge. Although this method is crude it should pro- 
vide an adequate representation of a model, where one 
assumes that BH merging is very efficient once the two 
BHs reside in the same parent galaxy. 

3. SIMULATIONS 

3.1. Galaxy models 

Our progenitor galaxies are setup using the method 
described by Springel et al. ( 2005b). The virial mass and 
radius of each model is determined by the virial velocity 
Uvir of the halo assuming a virial overdensity of Avir = 
200 using the following relations 



lOGHn ' 



lOGHo 



(1) 
(2) 



where Hq = 71 kms^^Mpc^^ is the present day Hub- 
ble parameter. Thus Wvir sets the mass of the model. 
The virial radius together with th e concentratio n param- 
eter c = Tvir/^s for a NFW halo (Nav arro et al.. 1997) is 
then used to construct iHerng uis9 ( 19901 ) pr ofile dark mat- 
ter halo s using the conversion detailed in I Springel et al.l 
(|2005bl ). For all galaxy models we assume a concentra- 
tion of c = 9, which is a typical value for a halo at 
z = (|Bullock et al.l[200l . All of our model galaxies 
thus correspond to loca l z = disk galaxies and unlike 
[Robertson et al.l (|2006cl ) we do not attempt a redshift- 
dependent scaling of our initial conditions. 

The dark matter halos are then populated with expo- 
nential disks with a constant mass fraction of rud = 0.041 
of the total virial mass resulting in a total disk mass of 
Md = mdM^iT- with a fractional gas content of /gas and 
the rest being st ars. The d isk sc ale length is deter- 
mined using the iMo et &\\ (|1998l ) formalism under the 
assumption that the fractional disk angular momentum 
jd equals the disk mass fraction md for a constant halo 
spin of A = 0.033 for all models. This assumption of 
Jd — jdJhaio {jd = md) correspouds to the conservation 
of the specific angular momentum of the material that 
forms the disk. The value for the spin parameter is close 
to the mean value f ound in numerical simulations (e.g. 
IVitvitska et al.ll2002i) . The vertical scaleheight zq of the 
stellar disk is radially constant and set to 0.2rd and the 
radial velocity dispersion is set to equal the vertical ve- 
locity dispersion. The vertical scale height of the gaseous 
disk is set via an integral constraint on the surface mass 
density, a self-consistently determined galaxy potential 
and the effective equation of state (EOS) of the mul- 
tipha se interstell ar medium model fMcKce fc Ostrikeil 
119771 Efstathio^ l2000t iJohansson fc Efstathioul l2006f). 
For the EOS we adopt the strongl y pressurized multi- 
phase model of the ISM, gEOS = 1 (.Springel fc Hernguisfj 
|2003|) that allows the construction of stable galaxy mod- 



els eve n with very large gaseous fractions (jSpringel et al.l 
l2005hD . 

All our galaxy models include a lHernguistl (|1990D stel- 
lar bulge component with constant total mass fraction 
of mi, = 0.01367 {Mb — mf,Mvh-), which is a third of 
the disk mass fraction. The bulge scale length b is also 
set to 0.2rd- Finally we insert a seed black hole with 
mass A/bh = lO^M© at rest in the center of each galaxy 
model. 

3.2. Simulation parameters 

At our standard numerical resolution we set the pri- 
mary galaxy models up with 20,000 gas and stellar disk 
particles, 10,000 bulge particles and 30,000 dark mat- 
ter particles. For the secondary galaxies in 3:1 merg- 
ers we use a third of the particle resolution of the pri- 
mary galaxies, thus resulting in 6,667 gas and stellar 
disk particles, 3,333 bulge particles and 10,000 dark mat- 
ter particles. We set the gravitational softening length 
of gas, newly formed stars and the black hole parti- 
cles to e = 0.1/i~^kpc and scale the softening lengths 
of the disk, bulge and more massive dark matter par- 
ticles with the square root of the particle mass result- 
ing in e = 0.2/i~^kpc for the bulge and disk particles 
and e = 0.8/i~^kpc for the dark matter particles. This 
ensures that the maximum gravitationa l force exerte d 
from a part i cle is in dependent of its mass (lDehnenll200l . 
iNaab et al.l (|20063) also showed that small variations in 
the respective softening lengths have only a minor effect 
on the central gas content, with typically only a 10% 
increase in the central gas content when equal soften- 
ing lengths and mass resolution were employed for all 
collisionless particles. The SPH properties of the gas 
and BH particle are averaged over the usual GADGET-2 
spline SPH kernel using ~ 64 SPH gas particles. Fur- 
thermore we require that the minimum SPH smooth- 
ing length is equal to the gravitational softening length. 
This choice safeguards against artificial stabilization o f 
small clumps at low resolution (|Bate fc Burkertjri997f) . 
In order to resolve the flow of gas around the black hole 
accurately we ran all simulations with high force accu- 
racy of Q!forco = 0.005 and tim e integration accuracy of 
7/acc = 0.02 (see lSpringeill2005l for details). 

The parameters governing the multi-phase feedback 
model are set as follows: The star formation timescale 
t'l — 8.4 Gyr, the cloud evaporation parameter Aq — 
4000 and the supernova 'temperature' Tsn = 4.0 x 10* K 
that reflects the heating rate from a population of su- 
pernovae for a given IMF. Thus, we here again follow 
ISpringel et al.l (|2005H ) who showed that these choices 
result in a star formation rate of ^ IMoyr^^ for a 
Milky Way-like galaxy. However, we note that these 
choices yield a comparatively low star formation rate 
that lies slightly below the observed Kennicutt relation 
(jKennicuttlll998f l. 

We set the dimensionless efficiency parameter of accre- 
tion in Eq. lAll to a — 100, which is considerably higher 
than the theoretically expected value of a ~ 1. This dis- 
crepancy is essentially due to limits imposed by the nu- 
merical resolution and the effective subresolution model 
for star formation. The estimated values for the gas den- 
sity and thermal sound speed are large-scale averages of 
the gas at the location of the BH resulting in artifically 
low densities and high sound speeds compared to the case 
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Fig. 1. — The star formation rate (top), the black hole accretion 
rate (middle) and the evolution of the black hole mass (bottom) 
as a function of time for four isolated galaxy models with initial 
gas mass fractions of 20% (black), 40% (green) and 80% (red) and 
virial velocities as indicated in the Figure. 

where the multi-phase structure of the ISM could be spa- 
tially resolved on the scale of the Bondi radius. Thus the 
larger value for a can essentially be seen as an empiri- 
cal correction factor that translates from the resolvable 
low mean density to the time-averaged small-scale den- 
sity at the location of the BH. A s imilar conclusion was 
reached bv iKhalatvan et all (|2008f) . who set a = 300 in 
their simulations. The exact choice of a within a fac- 
tor of a few is not critical as long as a is large enough 
to allow a black hole in a gas-rich environment to reach 
the Eddington regime, (again see iSprineel et al]l2005bl 
for details). As detailed in [JX]we set the radiative effi- 
ciency to = 0.1 and the thermal coupling constant of 
the black hole feedback energy to e/ = 0.05. 

All simulations presented in this paper were evolved 
for a total of t = 3 Gyr using the local Altix 3700 Bx2 
machine hosted at the University Observatory in Munich. 

3.3. Isolated disks with BHs 

We begin by running a representative sample of our 
disk galaxy models in isolation in order to study the sta- 
bility of the constructed galaxy models. The models are 
run at our standard numerical resolution. In Fig. [1] we 
show the resulting star formation rates, BH accretion 
rates and BH mass growth for four models with 20% 
(black), 40% (green), 80% (red) initial gas mass frac- 



TABLE 1 

Isolated galaxies and the A/bh — c and A/bh — A/, relations 



! ■ ^ 
vir 


/gas 


A/bh' 


•^bulgc 


A/bulgc'' 


^.^BH.cr 




80 


0.4 


5.12 


57.0 


0.23 


0.59 


0.22 


111 


0.2 


6.50 


78.3 


0.61 


0.21 


0.09 


160 


0.2 


56.0 


120.4 


1.83 


0.32 


0.24 


320 


0.8 


2640 


252.9 


14.7 


0.76 


1.08 



"Virial velocity in [kms ^] 

'Final mass of the BH in W^Mq 

"Final bulge stellar velocity dispersion in km/s 

''Final bulge mass in Mq 

"The ratio of the final BH mass to the BH mass expected from 
Eq. [3] if the final BH would lie on the A/bh ~ relation. 

^The ratio of the final BH mass to the BH mass expected from 
Eq. [4] if the final BH would lie on the A/bh — A/, relation 

tions. All galaxies show stable evolution in their star 
formation and BH accretion rates, reproducing the iso- 
lat ed disk results f or a sti ff EOS with (/eos = 1 presented 
bv lSpringel et all (|2005b[ ). 

All models start with a seed black hole with mass 
Mbh — IO^Mq at rest in the center, which grows due 
to gas accretion during the simulation as is shown in the 
bottom panel of Fig. [TJ After completing the isolated 
disk simulations we calculate the final black hole mass of 
the galaxies together with the mass- weighted line-of-sight 
stellar velocity dispersion Cbuigc measured from the bulge 
stellar particles using 50 randomly projected realizations 
of the galaxy models. Using the extracted velocity dis- 
persions we then calculate the expected BH masses of the 
galaxies and check if they lie on the observed Mbh — cr 
relation, 

log(A/BH/Af0) + 61og(CT/cro), (3) 

where the relation is defined relative to ( Tq = 200 kms~^. 
The o bserved values as determined by iTremaine et al.l 
l|200l are a — 8.13 ± 0.06 for the normalization coeffi- 
cient and h — 4.02 ± 0.32 for the slope. We also evaluate 
the total bulge stellar mass A/buigc by adding up the mass 
of all the stellar bulge particles. By definition this mass 
is not changing throughout the simulation, as the newly 
formed stars form from gas in the disk and we do not 
include their potential contribution to the final stellar 
bulge mass. We then compare the final BH masses with 
the expected BH masses if the galaxies would lie on the 
observed Mbh — M* relation, where we now use M* to 
denote the total bulge mass. 

log(MBH/M0) - c + dlog(Af,/10"MQ). (4) 

The best fi t ob served values as determined by 
iHaring fc R b? (2004) are c = 8.20 ± 0.10 for the nor- 
malization coefficient and d = 1.12 ± 0.06 for the slope. 
The final BH masses, stellar velocity dispersions, final 
stellar masses and the ratios between the final BH mass 
and the BH mass derived from Eqs. [3] and [5] assuming 
that the galaxies lie on the Mbh — c and A/bh ^ A/* 
relations respectively are tabulated in Table [TJ 

The initial seed mass of Mbh = IQ^Mq corresponds 
to a velocity dispersion of cr = 33 kms~^ (Eq. [3]), 
which would be the corresponding velocity dispersion of 
a galaxy with Wvir ^ 50 kms^^. Thus the BHs of the 
galaxies are not initially on the A/bh — c relation, but 
evolve towards this relation during the simulation. The 
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bulge velocity dispersions of the galaxies increase only 
slightly on the ^ 10% level during the simulations. The 
governing factor determining the evolution of the BHs 
for the isolated disk galaxies is the depth of the poten- 
tial well and the gas mass fraction. This can be seen in 
the ratio of the final BH mass to the expected BH mass 
from Eq. [3] tabulated in Table [T] The lower the gas mass 
fraction the lower the final BH masses are with respect 
to the expected Mbh — relation. Models with lower 
gas mass fractions would require additional gas inflow 
to the central regions in order to ensure sufficient BH 
growth. Such additional gas inflow could be triggered 
by galaxy mergers. The same initial BH seed mass of 
Mbh = IO^Mq would require initially a total bulge mass 
of = 1.4 • IO^Mq. Again this would only be vahd for 
the very lowest mass initial galaxy models and thus gen- 
erally the galaxies do not initially lie on the Mbh — M^. 
relation. Typically the simulations of the isolated galax- 
ies result in BH masses that are too low for their given 
stellar bulge masses. Again an additional central gas 
flow, potentially induced by a galaxy merger is required 
for the central BHs to reach the masses expected by their 
corresponding stellar bulge masses. 

Interestingly the simulation of the most massive and 
gas-rich galaxy (the Las — 80%, Wvir ^ 320 kms~^ model 
in Fig. [1] and Table results in a final BH mass that 
is relatively close to the BH mass expected from the 
-^BH — c relation and even superceeds the BH mass ex- 
pected from the Mbh — relation. However, some of 
this BH mass growth may be artificially high, due to 
the fact that the current BH feedback model does not 
restrict the size of the BH accretion radius. Thus the 
accretion radius of a a massive BH in a gas-rich region 
may grow to unphysically large sizes and accrete mate- 
rial that would physically not end up in the BH because 
of its high angular momentum and large distance from 
the BH. This is typically not a problem in merger simu- 
lations, where most of the BH mass grows by accretion 
of nearby high-density gas. 

3.4. The effect of resolution on the A/bh — cr relation 

We then proceed with a detailed analysis of the effect 
of the BH merger prescription and resolution on the fi- 
nal mass of BHs and the corresponding stellar velocity 
dispersions in 3:1 and 1:1 merger remnants. We simu- 
late 3:1 mergers with a quarter (tiny), half (low), double 
(medium) and four times (high) the number of particles 
compared to our standard numerical resolution ( §3.2|) . 
In addition we compute a control sample of 1:1 mergers 
at the standard, low and tiny resolutions. The five 3:1 
and three 1:1 initial conditions arc simulated using both 
the standard BH merging prescription and the BH repo- 
sitioning method described in i i2.2l resulting in a total 
of 16 simulations (see Table [2] for details). The gravita- 
tional softening lengths for all the particles including the 
BH are scaled with the inverse cube root of the number 
of particles, with larger numerical resolution resulting in 
smaller gravitational softening lengths. In addition we 
run a 3:1 and a 1:1 merger with a lower accretion effi- 
ciency of a = 1 at our standard resolution using the BH 
repositioning method. 

All galaxy models have a 20% initial gas fraction, with 
Vvir,p — 160 kms~^, Wvir,s = HI kms^^ and are setup 
with the parameters described in H3.1\ (p and s denote the 
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1:1 Repos 




^ • 1:1 a=l 




A 3:1 Std 




A 3:1 Repos 




A 3: 1 a= 1 












A • 


70 90 


120 150 180 220 



Dispersion [krri/s] 

Fig. 2. — The A/bh — f relation as a function of resolution 
and merger prescription, with larger/smaller symbols indicating 
higher/lower resolution, respectively (Table [2} ■ Triangles show 3:1 
and circles 1:1 merger remnants. Black symbols employ the stan- 
dard merger prescription and the red symbols use the BH repo- 
sitioning method. The blue symbols also employ the BH reposi- 
tioning method but with a reduced accretion efficiency of a = 1. 
The solid/open symbols indicate cases where the BHs did/did not 
merge. (Ta ble [2t. The lin es show the observed relation with errors 
by Trcmain eetall pOOl ). 

primary and secondary galaxies, respectively). In order 
to isolate the effects of resolution and merging prescrip- 
tion all 3:1 and 1:1 mergers are run with identical merg- 
ing geometries and initial disk orientations. For the orbit 
and initial orientation we adopt geome try G13 for both 
the 3: 1 mergers and 1 : 1 mergers from iNaab fc BurkertI 
(|2003l ). The G13 geometry corresponds to the inchna- 
tions ip = —109, is = 180 and the arguments of the peri- 
center = 60, ujs — for the primary and secondary 
galaxies, respectively. The galaxies approach each other 
on a parabolic orbit where the initial separation of the 
progenitors is i?init — 0.5(rvir,p + fvir,s) and the peri- 
centric distance is rpori — T'd.p + ''d.s, where r-vir.p, J^d.p 
and rvir,s7 ?'d,s are the virial and disk scale radii for the 
primary and secondary galaxies, respectively. We note 
that radial parabolic or near-parabolic orbits are gen- 
erally motivated by stati stics from N-body simulations 
(|Khochfar fc Burkertl[200Ql . 

After completing the merger simulations we calculate 
the final black hole mass together with mass-weighted 
linc-of-sight stellar velocity dispersion a measured from 
all stellar particles within the projected half-mass ra- 
dius To using 50 randomly projected realizations of the 
merger remnant. This determination of a is similar to the 
method used by observers who measure the luminosity- 
weighted velocity disp ersion within an effective radius 
(|Gebhardt et al.ll2000f ). The simulated final BH masses 
and stellar velocity dispersions together with the time of 
the BH merger are tabulated in Table [2l 

For the 3:1 mergers, the merging prescription has a 
strong effect on the final BH mass of the merger remnant, 
with the standard prescription producing BH masses that 
are a factor of 2-3 lower than the corresponding BH 
masses using the repositioning method. The final stel- 
lar velocity dispersion on the other hand converges with 
increasing resolution, with the repositioning method re- 
sulting in marginally lower velocity dispersions. The 3:1 
mergers with the standard BH merger prescription result 
in BH masses that are too low for their corresponding 
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TABLE 2 

3:1 AND 1:1 BINARY BH MERGER RESOLUTION AND MERGER PRESCRIPTION STUDY 



Resolution'^ 


Mass ratio 


be 

tgas 


^diak 












A/BH.fin" 






Tiny Std 


3:1 


0.16 


0.32 


1.28 


5000 


2500 


7500 


31 


46.0 


136.9 


2.02 


Tiny Rcpos 


3:1 


0.16 


0.32 


1.28 


5000 


2500 


7500 


31 


118 


132.2 


2.81 


Low Std 


3:1 


0.125 


0.25 


1.00 


10000 


5000 


15000 


15.5 


57.0 


134.4 


2.00 


Low Rcpos 


3:1 


0.125 


0.25 


1.00 


10000 


5000 


15000 


15.5 


118 


127.7 


2.50 


Standard Std 


3:1 


0.1 


0.2 


0.8 


20000 


10000 


30000 


7.75 


81.7 


134.7 


— 


Standard Repos 


3:1 


0.1 


0.2 


0.8 


20000 


10000 


30000 


7.75 


126 


126.8 


1.82 


Medium Std 


3:1 


0.08 


0.16 


0.64 


40000 


20000 


60000 


3.87 


26.7 


133.5 


1.84 


Medium Repos 


3:1 


0.08 


0.16 


0.64 


40000 


20000 


60000 


3.87 


96.9 


120.3 


1.72 


High Std 


3:1 


0.0625 


0.125 


0.5 


80000 


40000 


120000 


1.94 


26.1 


129.8 


2.81 


High Rcpos 


3:1 


0.0625 


0.125 


0.5 


80000 


40000 


120000 


1.94 


54.2 


129.5 


1.86 


Tiny Std 


1:1 


0.16 


0.32 


1.28 


5000 


2500 


7500 


31 


420 


155.4 


2.21 


Tiny Rcpos 


1:1 


0.16 


0.32 


1.28 


5000 


2500 


7500 


31 


508 


150.7 


1.57 


Low Std 


1:1 


0.125 


0.25 


1.00 


10000 


5000 


15000 


15.5 


134 


175.4 


1.84 


Low Rcpos 


1:1 


0.125 


0.25 


1.00 


10000 


5000 


15000 


15.5 


469 


153.1 


1.57 


Standard Std 


1:1 


0.1 


0.2 


0.8 


20000 


10000 


30000 


7.75 


456 


155.8 


1.72 


Standard Rcpos 


1:1 


0.1 


0.2 


0.8 


20000 


10000 


30000 


7.75 


516 


157.9 


1.50 



"•The secondary has a third of the particles of the primary 
^Gravitational softening lengths in /?,~^kpc 

tgas — Cstars — ^BH 
^disk ^bulgc 
■^Afgas = A^disk 

^Mass of a gas particle in 10^ Mq 

fFinal mass of the BH in lO^M© 

''Final stellar velocity dispersion within r^.^ in km/s 

'Time of BH merger in Gyr 

velocity dispersions as shown in Fig. [21 where the sim- 
ulated 3:1 (triangles) and 1:1 (circles) results are over- 
plotted on the observe d Mbh — cr relation from Eq. [3] by 
iTremaine et~all ()2002[ ). 

For the 1:1 mergers the effect of the merger prescription 
is less pronounced. With the exception of a single simula- 
tion at low resolution, all 1:1 mergers result in very simi- 
lar BH masses and velocity dispersions. At the standard 
resolution the two prescriptions produce final BH masses 
and velocity dispersions that agree within 10% with each 
other and the observed relation. All the 1:1 mergers with 
repositioning and two out of three of the standard pre- 
scription mergers lie on the observed Mbh — f relation 
as shown in Fig. [21 

Finally we show that reducing the accretion efficiency 
from our standard value of a = 100 to the physically 
motivated a — 1 produces BH masses that are far too 
low for their corresponding velocity dispersions (the blue 
points in Fig. [2). However, this is expected given the 
limitations imposed by the numerical resolution and sub- 
resolution modeling as was discussed in ij3.2l Changing 
the accretion efficiency or even removing the BH alto- 
gether does not affect the final stellar velocity dispersion, 
which is primarily set by the depth of the potential well. 
Thus, this implies that the BH is not dynamically im- 
portant for the galaxy as a whole. A similar problem 
would arise when comparing the simulated BH masses 
to the total stellar masses of the galaxies. Variations in 
numerical resolution and BH merging prescription affect 
only marginally the star formation rate and thus the BH 
masses in the standard prescription 3:1 merger remnants 
would be too low compared to their corresponding stellar 
masses. 

The numerical resolution tests show that the cen- 
tral stellar velocity dispersion converge within 10% with 
increasing resolution, whereas the BH mass converges 
within 25%, with the exception of the highest resolution 
simulation using the repositioning method. The highest 



resolution repositioning run produces a BH mass, which 
is lower by a factor of two compared to the other simu- 
lations. This is mainly due to the reduction of the gravi- 
tational softening of the BH particle. The accretion rate 
and mass growth of the BH is determined by a combi- 
nation of the accretion efficiency {a in Eq. lAip and the 
gravitational softening. Our chosen value of accretion 
efficiency was adopted for our standard numerical reso- 
lution, thus it is not surprising that smaller gravitational 
softening lengths produce slightly smaller BH masses. 

In Appendix [Cl we study in more detail the differences 
between the standard and repositioning merger prescrip- 
tion for both 3:1 and 1:1 binary galaxy mergers with 
BHs. We find that the repositioning merger prescription 
at the standard numerical resolution produces stable re- 
sults and thus we adopt this prescription and numerical 
resolution for all remaining simulations presented in this 
paper. 

4. RESULTS 

4.1. Unequal-mass 3:1 disk mergers with BHs 

Using the galaxy models described in §3.11 we set up 
five unequal-mass 3:1 mergers and five equal-mass 1:1 
mergers with three different initial gas mass fractions 
of /gas = 0.2, 0.4, 0.8, resulting in a total of thirty sim- 
ulations (the SI- and S2-series in Table [3]). The pro- 
genitor galaxies are at the standard numerical resolu- 
tion with the virial velocities of the models ranging from 
I'vir.p = 50 — 500 kms""'^ and Wvir.s = 35 — 347 kms~^ for 
the primary and secondary galaxies, respectively. The 
progenitors merge on prograde, nearly parabolic, in plane 
orbits (GO orbit), with inclinations ii = 12 = and argu- 
ments of the pericenter uji — uj2 — 0. The initial separa- 
tion and pericentric distances are set to the mean of the 
virial radius and twice the mean of the scale radius re- 
spectively, as detailed in §3.41 In addition we simulate a 
set of three 3:1 and three 1:1 mergers with varying initial 
disk geometries and orbits. These simulations (the Gl- 
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TABLE 3 

Unequal- and equal-mass mergers 



Models 




ir,proKl 






'ir,proK2 


Mass 






Orbit 






gas, tot 


"bulKo,tot 


^DM,tot 




Sl-l;l 


50, 80, 


160, 320, 500 


50, 


80, 


160, 320, 500 


1 


1 




GO 


0.2, 0.4, 


0.8 


40000 


20000 


60000 


15 


S2-3:l 


50, 80, 


160, 320, 500 


35, 


55, 


111, 222, 347 


3 


1 




GO 


0.2, 0.4, 


0.8 


26667 


13333 


40000 


15 


Gl-l:l 




160 






160 


1 


1 


G7. 


GIO. G13 


0.2 




40000 


20000 


60000 


3 


G2-3:l 




160 






111 


3 


1 


G4, 


GIO, G14 


0.2 




26667 


13333 


40000 


3 



"All virial velocities in kms 



AT 



disk 



1000.0 r 




1000.0 r 



1.0 1.5 3.0 

Time [Gyr] 



3.0 




1.0 1.5 3.0 

Time [Gyr] 



3.0 



Fig. 3. — The total star formation rate (top), the total black hole 
accretion rate (middle) and the evolution of the total black hole 
mass (bottom) as a function of time for three 3:1 (solid lines) and 
one 1:1 (dashed line) merger with initial gas mass fractions of 20% 
(black), 40% (green) and 80% (red). The filled circles indicate the 
time of merging of the BHs. 



Fig. 4. — The total star formation rate (top), the total black 
hole accretion rate (middle) and the evolution of the total black 
hole mass (bottom) as a function of time for three 3:1 mergers with 
varying initial gas mass fractions (black, green, red lines) and three 
3:1 mergers with varying orbital and initial geometries for a fixed 
gas fraction (blue lines). 



and G2-series in Table ^ were all run with /gas = 0.2 
and i;vir,p = 160 kms~^, Wvir,s = HI knis^^ progen- 
itors on the G7, GIO, G13 orbits for the equal-mass 
mergers and the G4, GIO, G14 orbits fo r the unequal- 
mass mergers fsee lNaab fc Burkertll2003l for a definition 
of the orbital parameters). The orbits in the G-series 
were chosen to provide a good coverage of the expected 
anis otropy-ellipticity (S — e) plane for merger remnants 
fsee lBurkert et al] l2007' for definitions of 6 and e). 

In Fig. [3]we plot the total star formation rates, BH ac- 
cretion rates and the BH mass growth for three 3:1 merg- 



ers (solid lines) with 20% (black), 40% (green), 80% (red) 
initial gas mass fractions together with one 20% gas frac- 
tion 1:1 merger (dashed lines). The star formation is very 
efficiently terminated by the BH feedback in 1:1 mergers, 
compared to a generally much shallower decline in star 
formation for 3:1 mergers. In the lowest mass 3:1 merger 
the BH feedback is actually unable to terminate the star 
formation and the galaxy is forming stars at a constant 
rate even after the BH merger. A similar trend is true 
for the BH accretion rates, with a much gentler decline 
in the accretion rate for galaxies undergoing 3:1 mergers 
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Fig. 5. — The Mbh — f (left) and Mbh — M, (right) relations for our complete sample of 3;1 (triangles) and 1;1 (circles) mergers. The 
black, green and red symbols show the effect of varying the initial gas mass fraction for our SI- (equal-mass) and S2-samples (unequal-mass). 
The blue symbols show the variation caus ed by the orbit and initial geometry fo r our Gl- and G2-samples (see also Table [Sj. The lines 
show the observed relations with errors bv lTremaine et al.l 11200^^ (left panel) and lHaring fc flTxi (2004) (right panel), respectively. 



TABLE 4 

Best fit Mbh — c relation for 3:1 and 1:1 mergers 



Sample N a b Aip^ 

Tot sample 36 8.07 ±0.06 3.82 ±0.15 0.29 

3:1 sample 18 8.06 ± 0.08 3.78 ±0.18 0.33 

1:1 sample 18 8.05 ± 0.07 3.77 ±0.18 0.26 

S1-S2 20% gas sample 10 7.85 ± 0.04 3.47 ±0.12 0.13 

S1-S2 40% gas sample 10 8.13 ±0.05 3.96 ±0.13 0.14 

S1-S2 80% gas sample 10 8.35 ±0.10 3.77 ±0.28 0.29 

Observed sample" 31 8.13 ±0.06 4.02 ± 0.32 0.25-0.3 

°Prom lTVemaine et al.] pOOl ). 

compared to 1:1 mergers. This evolution is mirrored in 
the growth of the BH masses, with the final BH mass be- 
ing lower in 3:1 mergers typically by a factor of 2-5, but 
with a relatively large scatter depending on the progen- 
itor masses and initial gas fractions. Most of the mass 
growth of the BHs typically takes place shortly after the 
merging of the BHs as indicated by the filled circles in 

Fig. m 

In Fig. |4] we show the results for three 3:1 mergers 
with varying initial gas mass fractions for a fixed orbit 
and initial disk geometry (black, green, red lines) and 
for three 3:1 mergers with varying initial orbits and ori- 
entations for a fixed gas mass fraction (blue lines). The 
initial gas fraction has a large effect on the height of the 
star formation and BH accretion peaks, with larger ini- 
tial gas fractions producing higher values, as expected. 
This results also in relatively large differences in the final 
BH masses, with the /gas = 0.8 simulations producing fi- 
nal BH masses that are larger by an order of magnitude 
compared to the /gas = 0.2 runs. The variation of the 
orbit and initial geometry for a fixed gas mass fraction 
produces much smaller differences. The peaks of the star 
formation rates and BH accretion rates only vary within 
a factor of two with changing orbits and initial disk ge- 
ometries. And although the time of the BH merging 
varies somewhat between different geometries the final 
BH mass is virtually unchanged. This was already shown 
to be the case for 1:1 mergers by'SpringeieO] (l2005bD 
and Fig. [4]confirms the same conclusion for 3:1 mergers. 
Hence the final BH mass is sensitive to the initial gas 
mass fraction of the galaxies, but relatively insensitive 
to the initial orbit and orientation of the galaxies. 



Figure [S] shows the Mbh — cr (left panel) and Mbh — 
M* (right panel) relations, extracted from our total 
sample of 36 simulated mergers, with the filled trian- 
gles and circles depicting 3:1 and 1:1 merger remnants, 
respectively, together with the observed relations by 
Tremaine et al. (2002) and Haring & Rix (2004). We 
perform a Levenberg-Marquardt least-squares fit to the 
Mbh — c relation as defined in Eq. [H Table [J] lists the 
best fit values for the normalization coefficient a and the 
slope b together with their la uncertainties derived using 
the Bootstrap method. In addition we list the dispersion 
^logMBH about each best-fit relation. The fitting anal- 
ysis was performed for our total merger sample, for the 
3:1 and 1:1 samples separately and for the sample split 
as a function of the initial gas mass fraction. Our best 
fit normalization coefficient and slopes for the total sam- 
ple are margin ally lower compared to the observational 
values of Tremaine et all (|2002t l. but still well within the 
error limits. Furthermore the intrinsic Mbh ~ cr disper- 
sion induced primarily by varying the initial gas frac- 
tion is in good agreement with the observed dispersion 
of AiogAfBH = 0-25 — 0.3. We also note that our results 
are in good ag reement with the z = Mbh — cr rela- 
tion derived bv [Robertson et all (|2006c| ). We do not find 
any statistically significant difference between the 3:1, 1:1 
and the complete merger sample, with all three samples 
resulting in virtually the same normalization coefficients 
and slopes for the Mbh — c relation. The normaliza- 
tion coefficient steadily increases with increasing initial 
gas mass fraction. The dispersion within a sample with 
constant gas mass fraction is in general low, with the 
exception of the sample with the largest gas mass frac- 
tion. This can be taken as a further indication that the 
dispersion in the Mbh — cr relation is primarily driven 
by changing initial gas mass fractions. We conclude that 
the simulated Mbh — c relation is overall robust and in 
agreement with the observed relation over a wide range 
of progenitor masses and initial gas fractions for both 
equal- and unequal-mass disk mergers. 

Similarly, we perform a least-squares fit to the A/bh — 
M* relation as defined in Eq. [H with the best fit param- 
eters, parameter errors and dispersions for our samples 
listed in Table \5\ We here define M* as the total mass 
of all luminous particles within the effective radius re, 
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including disk, bulge and newly formed stellar particles. 
This definition assumes that the merger remnant is a 
pure bulge system within this radius. We tested this 
assumption by performing a simple disk-bulg e decompo- 
sition following the methodology outlined in Abad i et al.l 
([2003) . We calculated the orbital circularity parameter 
ej = Jzl Jc\vc{E) for all stellar particles within the ef- 
fective radius, where Jz is the z-component of the spe- 
cific angular momentum and Jchc{E) is the specific an- 
gular momentum expected for a corotating circular orbit 
with binding energy E. We found that the fraction of 
particles on circular orbits, as defined by ej > 0.8 was 
well below 1% for 1:1 merger remnants and at most 2% 
for 3:1 merger remnants. Thus, for 1:1 and 3:1 mergers 
the assumption that the remnants are bulge-dominated is 
valid, whereas for higher mass ratio mergers this method 
might induce some additional scatter due to an underly- 
ing surviving disk comp onent in the final merger remnant 
(jNaab fc Truiillol 120061 ). In addition, we found that this 
definition for M* correlates reasonably well with total 
bulge mass as defined by the stellar bulge particles and 
finally thi s definition for is s i milar to the one em - 
ployed by iDi Matteo et all (|2008l ): ISijacki et all (|2007l ) 
and allows us to directly compare our results to theirs. 
Thus we adopt the simple definition of the total stellar 
mass within the effective radius as a proxy for the total 
bulge mass for all subsequent measurements of Af,. 

The best fit value of the normalization coefficient (c) 
for our total sample is in excellent agreement with the 
observations, but the derived slope d is too steep com- 
pared to the observed sample. The too large value of 
the slope is a general problem in all our simulation sam- 
ples as can be seen in Fig. [SJ For low-mass stellar 
systems the simulations typically produce too low BH 
masses, whereas for the higher-mass stellar systems the 
BH mass is slightly overpredicted, thus producing a tilt 
in the simulated relation with respect to the observed 
relation. The derived dispersion in the simulated rela- 
tion is marginally higher than the observed scatter, with 
the scatter being again primarily driven by variations 
in the initial gas mass fraction. Interestingly the final 
luminous mass within the effective radius remains vir- 
tually constant with changing initial gas mass fraction. 
In higher gas mass fraction simulations the lower initial 
stellar disk mass is compensated by higher star formation 
rates resulting in roughly the same final stellar mass. On 
the other hand, higher initial gas mass fractions produce 
larger final BH masses, thus resulting in vertical scat- 
ter in the Mbh — relation as can be seen in Fig. [5l 
This also results in a systematic increase in the normal- 
ization coefficient with increasing gas mass fraction, with 
the slope remaining relatively constant. A plausible ex- 
planation for the discrepancy between the simulated and 
observed slopes of the Mbh — -Af, relation lies in the em- 
ployed stellar feedback model. A more aggressive stellar 
feedback including galactic winds (|Springel &: Hernguistj 
l2003t lOpp enhciraer & Dave' ^GOG") would primarily affect 
the low-mass stellar systems, reducing their final stellar 
mass, thus lowering the simulated slope and bringing it 
in better agreement with the observed slope. I ndeed , 
both the cosmologi cal studies bv lDi Matteo et all (|2008f ): 
ISiiacki et al.l ()2007f ) found a better fit to the expected re- 
lations at the low mass end when they included a model 
for supernova-driven galactic winds in their simulations. 



TABLE 5 

Best fit Mbh — ^I* relation for 3:1 and 1:1 mergers 



Sample 


N 




c 






d 






Tot sample 


36 


8.17 


± 


0.10 


1.40 


± 


0.07 


0.44 


3:1 sample 


18 


8.04 


± 


0.11 


1.34 


± 


0.08 


0.47 


1:1 sample 


18 


8.24 


± 


0.10 


1.41 


± 


0.10 


0.38 


S1-S2 20% gas sample 


10 


7.86 


± 


0.07 


1.34 


± 


0.05 


0.17 


S1-S2 40% gas sample 


10 


8.28 


± 


0.08 


1.45 


± 


0.06 


0.22 


S1-S2 80% gas sample 


10 


8.68 


± 


0.13 


1.36 


± 


0.12 


0.29 


Observed sample" 


30 


8.20 


± 


0.10 


1.12 


± 


0.06 


0.30 



"From lHarine fc Rixl 1 12001 ). 



An alternative explanation for the discrepancy lies in the 
fact that the progenitors of lower mass galaxy mergers 
could have larger initial gas mass fractions compared to 
the higher mass merger progenitors. This would also ef- 
fectively remove the tilt in the derived correlation bring- 
ing it in better agreement with the observed relation. 
Given these limitations in the stellar feedback model em- 
ployed in our simulations we find that the final fit to the 
Mbh — relation is adequate, with a well defined nor- 
malization and a slightly steeper slope compared to the 
observational relation. 

4.2. The accretion history of BHs in unequal mergers 

We now proceed to study the BH accretion and star 
formation histories for unequal-mass mergers with vary- 
ing mass ratios. To this end we ran five mergers with 
mass ratios of 1:1, 2:1, 3:1, 4:1 and 6:1 on co-planar pro- 
grade orbits (GO). This orbital configuration was cho- 
sen in order to maximize the non-axisymmetric torques 
during the merger resulting in the most inten sive bursts 
in BH accretion and star formation activity (|Cox et al.l 
l2008l ). All the progenitor galaxies had initial gas mass 
fractions of 20% with the primary progenitor having 
Wvir = 160 kms~^ at the standard numerical resolution. 
The secondary galaxies were scaled down in mass and nu- 
merical resolution accordingly in order to maintain equal 
mass resolution per particle. The initial seed BH mass 
was set to MBH,init = lO^M© for all models. Again, the 
initial separation and pericentric distances are set to the 
mean of the virial radius and twice the mean of the scale 
radius respectively as detailed in H3A[ 

In Fig. [6] we show the evolution of the resulting star for- 
mation (top), BH accretion (middle) and BH mass (bot- 
tom) for the five mergers as a function of time. Increasing 
the mass ratio of the merger systematically lowers the 
peak star formation rate and increases the duration of 
star formation activity after the merger. For the highest 
mass ratio merger the star formation rate is virtually con- 
stant throughout the simulation with only a mild peak 
during the first passag e (see also idi Matteo et al.l 120071 : 
iBournaud et al.|[2007bl ). A similar evolution is seen in the 
BH accretion rates, with the BH accretion rates system- 
atically peaking at higher rates with decreasing merger 
mass ratios. For the higher mass ratio mergers no well 
defined BH accretion peak is discernable, instead the BH 
accretion rates stay roughly constant on an elevated level 
for a prolonged time period compared to the pre-merger 
BH accretion rates. The evolution of the BH accretion 
rates are mirrored in the growth of the BH masses, with 
the final BH mass being systematically lower for increas- 
ing mass ratio of the merger. Furthermore, the slope 
of the A/bh growth as a function of time becomes shal- 
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Fig. 6. — The total star formation rate (top), the total black hole 
accretion rate (middle) and the evolution of the total black hole 
mass (bottom) as a function of time for co-planar prograde (GO) 1:1 
(black), 2:1 (green), 3:1 (red), 4:1 (blue) and 6:1 (orange) mergers. 
The filled circles indicate the time of the BH merger, the filled 
diamonds in the top panel and the filled triangles in the bottom 
two panels show the location of the maximum star formation and 
BH accretion peaks, respectively. 

lower with increasing merger mass ratio. There is also 
a systematic delay in the time of the BH merger with 
increasing mass ratio, as indicated by the filled circles in 
Fig. [6l For the lower mass ratio mergers the peak of 
the BH activity (the filled triangles in Fig. (S]) typically 
occurs shortly after the merging of the BHs, whereas for 
the higher mass ratio mergers the peak of the BH activity 
is not directly related to the merging time of the BHs. 

In Fig. [7] we study the duration of the BH accretion 
and star formation activity as a function of merger mass 
ratio. We define a half-mass growth time AT1/2 during 
which half of the final BH mass and half of the total 
stellar mass is formed respectively. In both cases Ari/2 
is centered on the peak of the corresponding activity, 
the BH accretion on the maximum of the BH accretion 
rate (the triangles in Fig. (5]) and the star formation rate 
on the peak of the SFR marked with diamonds in Fig. 
[6l The resulting half-mass timescales are shown in the 
bottom panel of Fig. [71 The ATi/2gjj is strongly cor- 
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Fig. 7. — The maximum BH accretion and star formation rates 
and timescales as a function of merger mass ratio. The bottom 
panel shows the half mass growth time of the BHs (circles) and the 
stellar mass (triangles) centered on the maximum BH accretion 
and star formation rates, respectively. The top panel shows the 
corresponding peak BH accretion (circles) and star formation rates 
(triangles) during the corresponding half mass growth time 
The dashed line is the best resulting fit of Eq. [5]to the peak values. 



related with the mass ratio of the merger. For 1:1 and 
2:1 mergers the growth of the M-qu is very concentrated 
in time, with half of the final BH mass growth occur- 
ring in less than 100 Myr. For the 3:1 and 4:1 mergers 
ATi/2bjj ^ 0.5 Gyr, with the 6:1 merger resulting in 
ATx/2gjj ^ 1 Gyr. In the lower mass ratio mergers a 
significant fraction of the gas is funnelled to the center 
during the merger, enabling the BHs to experience rapid 
Eddington limited growth (Eq. IA2p . For larger mass ra- 
tio mergers the tidally driven non-axisymmetric torques 
are weaker and a lower fraction of the gas reaches the 
center resulting in lower BH accretion rates and longer 
BH half-mass growth timescales. 

The corresponding stellar half-mass timescales 
ATi/2j,jg^j. (triangles in the bottom panel Fig. [7]) also 
show a clear correlation with the mass ratio of the 
merger. In the 1:1 merger half of the final stellar 
mass is formed in a short major burst lasting about 
ATi/2 , . ~ 0.15 Gyr. This value is comparable to 
the s tar formation timcscale derived by ICox et al.l 
l)2008f ). who calculated a full width at half maximum of 
FWHM ~ 0.1 Gyr for the star formation peak of a typical 
1:1 merger. For higher mass-ratio mergers the resulting 
star formation timescales are ~ 0.7 — 1.0 Gyr, 

with the highest mass-ratio merger having the longest 



timescale of AT1/2 



1.5 Gyr. In these cases it is 
not straightforward to relate our ATi/2 ^).^^ to the burst 
timescales defined by ' Cox et all (|2008l ). as the higher 
mass ratio mergers have less well defined star formation 
peaks with the star formation distributed more evenly 
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over time. However, the general trend that increasing 
the merger mass ratio results in increasingly longer 
star formation ti mescales agrees well with the results of 
iCox et all (I200I . 

Finally we plot in the top panel of Fig. [7] the cor- 
responding maximum black hole accretion and star for- 
mation rates during their respective half mass growth 
timescales. By defining the variable q as the mass ra- 
tio between the primary and secondary component we 
can fit the logarithms of the peak BH accretion and star 
formation rates with the following linear relation 

log Maximum [MQyr""'^] = ap -l- ai x q, (5) 

where aq and ai are the inferred normalization and slope 
respectively. Both the maximum BH accretion rates and 
star formation rates are well fitted by Eq. [5] (dotted lines 
in Fig. [7]) resulting in the following best fitting param- 
eters respectively: (ao^BHAR = 0-S^,ai^BHAB. = -0.47) 
{ao,SFR = 1-82, ai^sFfl; = —0.24). The ratio of the peak 
BH accretion rate to the peak star formation rate is thus 
of the order of MeH.peak/^SF^pcak ^ 10~^. On the other 
hand, the ratio of the mean BH accretion rate to the 
mean star formation averaged over the whole simulation 
is closer to MBH,moan/^SF,mean ~ 10~^, with this ratio 
varying systematically between MBH.moan/^SF.mean = 
2-10-3 (1:1 merger) and MBH,mcan/MsF,mca„ = 0.2-10-3 
(6:1 merge r). Interestin gly, recent observations of 2 ~ 2 
galaxies bv lDaddi et al.l (2007b,. ,a) also find indications of 
an universal ratio of Mbh/-^sf ^ 10^^ and this is also 
expected from the observed Mbh — M^, relation. Natu- 
rally the normalizations and to some extent the slopes of 
the derived BH and SFR activity relations are dependent 
on the initial gas mass fraction and orbital geometry of 
the system. However, the derived ratios between the BH 
accretion and SF rates should be valid and we find in our 
simulations that for any given orbital configuration and 
initial gas mass fraction the ratio of the peak BH accre- 
tion and SF rates is ~ 10"^, whereas the corresponding 
mean value is of the order of ~ lO^^. 

4.3. Mixed E-Sp mergers with BHs 

We now turn our attention to mixed E-Sp mergers, in 
which a gas-poor early-type galaxy merges with a gas- 
rich spiral disk galaxy. To this end we simulated to our 
knowledge for the first time a large sample of mixed 
E-Sp mergers. For early-type progenitors we use the 
merger remnants of 1:1 and 3:1 mergers discussed in ij4.11 
whereas the progenitor disk galaxies are setup using the 
disk models discussed in §3.11 Our mixed merger sample 
consisted of five 1:1 merger remnants that are merged 
with the primary disks of the original simulations (Sl-E- 
Sp) and another five runs in which 3:1 merger remnants 
are merged with the corresponding secondary disks of the 
original simulations (S2-E-Sp). This sample of ten mixed 
mergers were run on the same GO orbit as the original 
mergers. In addition, we simulated six E-Sp mergers of 
1:1 and 3:1 merger remnants with primary and secondary 
disks respectively (Gl-E-Sp and G2-E-Sp) with varying 
initial orientations and orbital geometries as indicated in 
Table [6l This results in a total sample of 16 mixed E-Sp 
mergers. The initial gas fraction was set to 20% in both 
the original simulations that produced the E-progenitors 




: Sl-E-Sp: l:l-P:Sp v^=3S0 I 

10*' r • n 



Sl-E-Sp: l:l-P:Sp T^=ieO - 

o 10" ■ • ^ 

^ \ S3-E-Sp: 3:l-S:Sp v^= 1 1 1 E 
M " , A - 

m 

^ lO'' r Sl-E-Sp: 1-1-P'3p v^=60 H 



10'' r 1 

10^ r . , , . I . . , , I , . . . I , , , , I . . . , I , , , , " 

0.0 0.5 1.0 1.5 3.0 2.5 3,0 
Time [Gyr] 

Fig. 8. — The total star formation rate (top), the total black hole 
accretion rate (middle) and the evolution of the total black hole 
mass (bottom) as a function of time for four E-Sp mixed mergers. 
The solid lines show mixed E-Sp mergers of 1:1 merger remnants 
with primary spiral disks and the dashed line shows a mixed E-Sp 
merger between a 3:1 merger remnant and a secondary spiral disk. 
The gas disks have initially /gas = 20% and the virial velocities 
as indicated in the Figure. The filled circles indicate the time of 
merging of the BHs. 

and in the disk progenitors. The black hole mass of the 
early-type progenitors is set by the result in the orig- 
inal simulation and is typically A/bh ^ 10® — IO^Mq, 
whereas the initial seed BH mass of the disk galaxies is 
always fixed at Mbh — 1O^M0. 

We plot the corresponding total star formation rates 
(top), the total BH accretion rates (middle) and the evo- 
lution of the total BH mass (bottom) in Fig. [8] for a 
sample of four mixed E-Sp mergers. The three solid lines 
show mergers between 1 : 1 merger remnants and primary 
disks and the dashed line a merger between a 3:1 merger 
remnant and a secondary disk, with the color scaling in- 
dicating the virial velocity of the disk progenitors. The 
combined star formation rate of the mixed E-Sp mergers 
is lower than in Sp-Sp (Fig. ^ mergers, due to the lower 
amount of cold gas available for star formation. The disk 
progenitor contains /gas = 20% of gas initially, whereas 
the early-type progenitors typically have an initial gas 
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TABLE 6 
Mixed E-Sp mergers 



Models Progenitor 1 '»vir,prog2° Mass ratio Orbit /gaa.init ^Vdjak^tot -^bulgc.tot A^DM.tot N^j, 

51- E-Sp Sl-l:l 50, 80, 160, 320, 500 2:1 GO 0.2 60000 30000 90000 5 

52- E-Sp S2-3:l 35, 55, 111, 222, 347 4:1 GO 0.2 33334 16666 50000 5 
Gl-E-Sp Gl-l:l 160 2:1 G7, GIO, G13 0.2 60000 30000 90000 3 
G2-E-Sp G2-3:l 111 4^1 G4, GIO, G14 0.2 33334 16666 50000 3 

"All virial velocities in kms~^ 




Fig. 9. — The Mbh — (left) and Mbh — Af» (right) relations for our sample of mixed E-Sp mergers. The merger remnants of the 
original 1:1 and 3:1 mergers are shown by black and blue filled circles and triangles for the S and G merger series, respectively (Table [Sj. 
The mixed merging of 1:1 merger remnants with primary spiral disks with /gas = 20% (l:l-P:Sp) are shown with filled red (Sl-series) and 
purple (Gl-series) circles. The mixed merging of 3:1 merger remnants with secondary spiral disks with /gas = 20% (3:l-S:Sp) are s hown 
with filled red (S2-series) and pu rple (G2-series) triangles. The lines show the observed relations with errors bv lTremaine et al.l l(200^ (left 
panel) and IHarine &: Rixl II2004I) (right panel) , respectively. 



TABLE 7 

Best fit Mbh — f relation for E-E and E-Sp mergers 



Sample N a b Aipg j./^ 

Progenitor sample 16 7.83 ± 0.04 3.53 ±0.11 0.16 
E-Sp Mixed sample 16 8.03 ± 0.04 3.55 ±0.12 0.13 
E-E Remerger sample 16 8.13 ±0.05 3.41 ±0.10 0.18 



ing on the bulge mass of the progenitor galaxy. Thus the 
BHs of the disk galaxies are not initially on the Mbh — c 
relation and they are also not plotted in Fig. [9] Using a 
least-squares fit and Bootstrap technique we derive the 
best fitting parameters for the Mbh — cr relation of our 
mixed merger and progenitor samples (see Table [7]). For 
the mixed mergers the normalization parameter a is 0.2 
dex higher than for the Sp-Sp progenitor sample, but 0.1 
dex lower than for the E-E mergers. The combined in- 
situ growth due to gas accretion and merging of the BHs 
thus result in final BH masses that arc very similar to the 
ones derived in both the Sp-Sp (Fig. [5]) and E-E mergers 
(Fig. [TT|) . The mixed merging of E-Sp galaxies seems to 
maintain the slope b of the Mbh — c relation compared 
to the progenitor sample and even slightly reduces the 
dispersion. The mixed mergers primarily increase the 
final BH mass, with smaller changes in the final veloc- 
ity dispersion compared to the initial dispersion of the 
early-type progenitors. Ov erall the slope a.nd dis persion 
are lower than the observed lTremaine et al.l (120021 ) values 
of 6 = 4.02 ± 0.32, Aiog = 0.25 - 0.3. This is primar- 
ily due to the fact that all of the original mergers and 
disk progenitors had initial gas mass fractions of 20%. 
Including mixed mergers of merger remnants and disks 
with higher initial gas mass fractions would increase the 
normalization, slope and scatter in the derived Mbh — cr 
relation (see Table |4]). The mixed mergers undergo a 
similar phase of strong self-regulated BH growth as the 
Sp-Sp mergers, allowing the BH feedback to regulate the 
final BH mass in agreement of the Mbh — cr relation. 
Thus we conclude that the Mbh — o" relation should re- 
main relatively stable even after several generations of 
mixed E-Sp mergers. 



mass fraction of /gas ^ 5%, with typically only ~ 1% of 
this gas being cold and dense. After the merger of the 
E-Sp galaxies the star formation rate declines rapidly in 
all the merger remnants. The black hole accretion rates 
are somewhat lower than in Sp-Sp mergers, but still rela- 
tively high due to the large amount of available gas in the 
central regions of the merger remnants originating from 
the disk progenitors. The BH accretion rates show a clear 
peak, both during the time of the first passage and at the 
time of final coalescence. We define a BH mass growth 
factor as /BH,insitu = AMBH,insitu/MBH,finai, which gives 
the ratio of BH mass growth due to gas accretion dur- 
ing the simulation with respect to the final BH mass. 
Quantitatively, the fraction of the BH mass that accu- 
mulates by gas accretion during the mixed mergers is in 
the range of /BH,insitu ~ 20 — 50%, with typical mean val- 
ues of /BH,insitu ~ 30%. Thus in most cases the majority 
of the final BH mass is still contributed by the massive 
seed BHs residing in the early-type progenitors. 

The Mbh - cr (left panel) and Mbh - M, (right panel) 
relations for our mixed E-Sp merger (red and purple sym- 
bols) sample together with the E-type progenitor popu- 
lation (black and blue symbols) is plotted in Fig. [9l All 
the progenitor disk galaxies have initial seed BH masses 
of A/bh = lO^M© with the velocity dispersion a depend- 
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TABLE 8 

Best fit Mbh — relation for E-Sp and E-E mergers 



Sample 


N 


c 


d 


^loR Mrh 


Progenitor sample 


16 


7.78 ±0.07 


1.35 ± 0.05 


0.21 


E-Sp Mixed sample 


16 


7.83 ± 0.05 


1.39 ±0.07 


0.16 


E-E Remerger sample 


16 


7.86 ± 0.05 


1.38 ± 0.04 


0.17 



Finally, we also perform a least-squares fit to the 
Mbh — Af* (Eq- SI) relation for our mixed E-Sp sample 
with the best fit parameters listed in Table [51 The mixed 
merging of an E-Sp galaxy population effectively main- 
tains the original A/bh — relation of the progenitor 
sample, even slightly reducing the intrinsic dispersion in 
the relation. The self-regulated stellar and BH feedback 
models regulate the star formation process and BH mass 
growth in mixed mergers similarly to Sp-Sp mergers thus 
ensuring that the derived original A/bh — relation is 
maintained. We thus conclude that even repeated gener- 
ations of mixed E-Sp mixed mergers should not greatly 
affect the observed ATbh — relation. 

4.4. E-E remergers with BHs 

Finally, we studied E-E mergers including BH feedback 
by remerging a sample of our merger remnants produced 
by the 3:1 and 1:1 mergers discussed in ij4.1l We set up 
and ran five equal-mass remergers of our Sl-l:l (Table 
[3]) merger remnants (Sl-E-E), with the initial progenitor 
disks of the original simulations having virial velocities 
in the range Uvir = 50 — 500 kms~^. In addition, we 
remerged five of our 1:1 merger remnants with the corre- 
sponding set of 3:1 merger remnants (S2-E-E), resulting 
in a mass ratio of 3:2 (see Table IH]). All of the mergers 
were run on identical orbits to the original mergers and 
the initial gas fraction of the progenitors for the orig- 
inal simulations was fixed at 20%. In addition to the 
ten S-series remergers that were run on the GO orbit, 
we simulated three E-E mergers of 1:1 merger remnants 
(Gl-E-E) and three E-E mergers of 3:1 merger remnants 
(G2-E-E) with varying initial orientations and orbital ge- 
ometries as indicated in Table [51 thus resulting in a total 
sample of 16 E-E remergers. 

In Fig. [TOlwe plot the star formation rate (top), the BH 
accretion rate (middle) and the evolution of the BH mass 
(bottom) for three 1:1-3:1 E-E remergers (solid lines) and 
one 1:1-1:1 E-E remerger (dashed line). Here the color- 
ing indicates the virial velocities of the primary galaxies 
in the original simulations, Uvir,p — 80 kms~^ (green), 
Vvir,p — 160 kms^^ (black) and Uvir.p = 320 kms~^ 
(red). The initial star formation rates are generally very 



low due to the low gas fractions of /g. 



1 - 5% and 



/gas ~ 5 — 10% for the 1:1 and 3:1 merger remnant E- 
progenitors, respectively. The initial gas fraction directly 
depends on the strength of the initial interaction that 
gave rise to the merger remnants used as progenitors for 
the E-E remergers. Increasing the masses of the progen- 
itors and using more direct planar orbits (GO) produces 
more violent initial encounters, thus decreasing /gas- In 
addition, most of the gas is hot with low densities and the 
fraction of cold star- forming gas is typically below < 1% 
of the total gas mass. As can be seen in Fig. [TUlthe star 
formation is very effectively terminated shortly after the 
merging of the E-E progenitors on comparable timescales 
to 1:1 Sp-Sp mergers and thus more efficiently than in 3:1 



10.0 r 




m 

10"" 
lO*' 



10' r 



10*' r 



10° 




Sl-E-E: 1:1-1:1 v. ,= 160 



SS-E-E: 1:1-3:1 v^j,= 160 
SS-E-E: 1-1-3:1 v„ =eO 



0.0 0.5 1.0 1.5 3.0 
Time [Gyr] 



2.5 



3,0 



Fig. 10. — The total star formation rate (top), the black hole 
accretion rate (middle) and the evolution of the total black hole 
mass (bottom) as a function of time for four E-E remergers. The 
solid lines show remergers between 1:1 and 3:1 merger remnants 
and the dashed line shows a E-E remerger between two 1 : 1 merger 
remnants. The virial velocities of the primary galaxies in the orig- 
inal simulations are as indicated in the Figure. The filled circles 
indicate the time of merging of the BHs. 

Sp-Sp and mixed E-Sp mergers. For the 1:1-1:1 E-E re- 
merger the star formation rate is zero throughout the 
simulation as the star formation was already completely 
terminated in the original 1:1 merger that served as a 
progenitor galaxy for the E-E merger (see Figs. [3l and 

The black hole accretion rates strongly correlate with 
the available amount of gas in the central regions and 
are thus typically one to two orders of magnitude lower 
for dry E-E mergers compared to the mixed E-Sp and 
Sp-Sp mergers. This is also reflected in the total BH 
mass growth plots (Fig. [10]) that remain virtually flat 
throughout the simulations. The growth of the final BH 
mass is primarily due to the merging of the BHs with 
the in-situ growth of the BH mass through gas accre- 
tion being relatively low. The exact fraction of BH mass 
grown in-situ (/BH.insitu) due to gas accretion during the 
simulation with respect to the final BH mass, depends 
on the amount of available gas initially. For remerging 
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TABLE 9 

E-E REMERGERS 



Models Progenitor 1 Progenitor 2 Mass ratio Orbit /gaB,init -^diBk,tot -^bulgc,tot -^DM.tot ^sii 

51- E-E Sl-l:l Sl-l:l TIT GO 02 80000 40000 120000 s" 

52- E-E Sl-l:l S2-3:l 3:2 GO 0.2 66667 33333 100000 5 
Gl-E-E Gl-l:l Gl-l:l 1:1 G7, GIO, G13 0.2 80000 40000 120000 3 
G2-E-E G2-3:l G2-3:l Ijl G4, GIO, G14 0.2 26666 53334 80000 3 




Fig. 11. — The Mbh — c (left) and Mbh — ^* (right) relations for our sample of E-E remergers. The merger remnants of the original 1:1 
and 3:1 mergers are shown by black and blue filled circles and triangles for the S and G merger series, respectively (Table|3ll- The remerging 
of 1:1-1:1 merger remnants are shown with filled red (Sl-series) and purple (Gl-series) circles, the remerging of 1:1-3:1 merger remnants 
(S2-series) are shown with red diamond s ymbols and final l y the remerging of 3:1 -3:1 remnants ( G2-se ries) by purple triangles. The lines 
show the observed relations with errors by ITremaine et al.l 1120021) (left panel) and IHaring fc Rixl II2004I ) (right panel) , respectively. 



of 1:1 merger remnants this fraction is typically very low 
at /BH,insitu < 1%, whereas for the remerging of 1:1-3:1 
merger remnants the fraction is higher at /eH.insitu ^ 
10%, with the highest fraction of /BH,insitu ~ 20% pro- 
duced for remerging of 3:1-3:1 merger remnants, as ex- 
pected. There is relatively large scatter in the /sH.insitu 
fraction as a function of the orbital configuration and 
the initial masses of the progenitor galaxies primarily 
set by the corresponding scatter in the initial gas mass 
fractions of the progenitor galaxies. All the fractions 
above are estimates for only first generation remergers, 
with the /BH,insitu ratio rapidly decreasing with increas- 
ing generations of remergers. The same is true for the 
star formation, E-E remerging effectively quenches the 
star formation (Fig. [T0|) and already a second genera- 
tion of dry merging is expected to show very weak, if 
any, signatures of residual star formation. 

In Fig. [11] we plot the Mbh — cr (left panel) and the 
Mbh — A/» (right panel) relations for our E-E remerger 
(red and purple symbols) sample together with the pro- 
genitor population (black and blue symbols). The best 
fitted Mbh — f parameters together with the Icr uncer- 
tainties for our E-E remerger sample are listed in Table[71 
The normalization parameter a is 0.3 dex higher for the 
E-E remerger sample compared to the progenitor sample 
and 0.1 dex higher compared to the mixed E-Sp sample. 
However the slope 6 is marginally shallower compared to 
both the progenitor and mixed merger samples. The E- 
E remerging of the merger remnants increases the final 
BH mass, but leaves the velocity dispersion more or less 
unchanged, resulting in an increase of a. The relative 
increase in Mbh is slightly larger for low-mass systems 
compared to higher-mass systems, resulting in a slightly 
lower slope h. This is especially true for the remergers 



on the GO orbit, whereas remergers on other orbits (the 
purple G-series in Fig. [TT|) also increase slightly the ve- 
locit y dispersion of the resu lting merger remnants (see 
also (Robertson et al.ll2006bl ). The overall lower value of 
the slope h is again due to the fact that only E-E merg- 
ers of progenitors with initial gas fractions of /gas = 20% 
in the original simulations were included in construct- 
ing the E-E sample. We thus conclude that dry E-E re- 
merging typically increases the normalization parameter 
a and lowers the slope 6 slightly, but a single generation 
of dry remerging is not able to destroy the Mbh — cr re- 
lation initially generate d by gas-rich mergers (see also 
[Robertson etalTl2006bf) . E-E remerging also increases 
slightly the dispersion in the fitted relation, but the vari- 
ation in the initial gas fraction still remains the dominant 
source of scatter. Unlike Sp-Sp and mixed E-Sp mergers 
the gas-poor E-E mergers only experience a weak self- 
regulatory BH growth phase, with this phase missing al- 
together for the gas-poorest E-E mergers. This lack of 
the self-regulating feedback phase can thus explain some 
of the larger scatter in the final E-E remerger Mbh — <^ 
relation. 

We also perform a least-squares fit to the Mbh — M* 
(Eq. |4]) relation for both the E-E remerger sample and 
the progenitor sample with the best fit parameters, pa- 
rameter errors and dispersions listed in Table [51 Dry 
E-E remerging effectively maintains the normalization 
coefficient (c) and slope (d) of the progenitor sample, 
with a very minor increase of both parameters that lie 
within the \a uncertainties of the corresponding pro- 
genitor parameters. Furthermore dry remerging roughly 
maintains the dispersion derived for the progenitor pop- 
ulation. This is expected since the in-situ growth from 
gas of both the BH and stellar mass are very low dur- 
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Fig. 12. — The binding energy correlation, which is statistically equivalent to the BHFP. In the left panel we plot the correlation for our 
complete sample of 3:1 (triangles) and 1:1 (circles) mergers, with the colors indicating variations in the initial gas mass fractions (black, 
green, red) and in the orbital geometries (blue). In the right panel we plot the corresponding relation for our mixed E-Sp (circles and 
triangles) and the E-E remerge r sample (squares, diam onds and downward pointing triangles). The lines show the observed relation with 
scatter from the compilation of lHopkins et ahl l)2007al) 



TABLE 10 

Best fit binding energy relation for 3:1, 1:1, E-Sp and E-E 

MERGERS 



Sample 


N 




e 






f 






3:1 sample 


18 


8.05 


± 


0.10 


0.79 


± 


0.04 


0.41 


1:1 sample 


18 


8.15 


± 


0.08 


0.81 


± 


0.05 


0.32 


S1-S2 20% gas sample 


10 


7.85 


± 


0.05 


0.76 


± 


0.03 


0.14 


S1-S2 40% gas sample 


10 


8.22 


± 


0.06 


0.84 


± 


0.03 


0.17 


S1-S2 80% gas sample 


10 


8.54 


± 


0.12 


0.79 


± 


0.06 


0.28 


E-Sp Mixed sample 


16 


7.91 


± 


0.06 


0.77 


± 


0.03 


0.15 


E-E Remerger sample 


16 


7.98 


± 


0.05 


0.76 


± 


0.02 


0.18 


Total sample 


68 


8.04 


± 


0.05 


0.79 


± 


0.02 


0.31 


Observed sample" 


38 


8.23 


± 


0.06 


0.71 


± 


0.06 


0.25 



"Prom lHopkins et al.l l|2007al ') 



ing the simulation. For the hmiting case of no available 
gas and no in-situ growth of the BH and stellar mass, 
the original Mbh — M^, relation would be perfectly con- 
served. Thus we expect that increasing generations of 
E-E mergers would not substantially increase the scatter 
in the observed Mbh — M^, relation. The discrepancy 
between the derived Mbh — M, parameters for the E-E 
remergers and the observed parameters would be reme- 
died by remerging merger remnants with higher initial 
gas mass fractions (increases c) and by employing more 
effective stellar feedback models (lowers d) as was already 
discussed in 

4.5. The Black Hole fundamental plane 

Recent studies bv lHopkins et al.l ()2007al lbl) have shown 
that the well established Mbh — cr and Mbh — M» rela- 
tions are only two projections of a more fundamental un- 
derlying relation. This so called Black Hole Fundamental 
plane (BHFP) is defined as Mbh o-3.o±o.3^o.43±o.i9 
or Mbh ~ ]\,.jo.54±o.n^2.2±o.5 .^^^ jg analogous to 
the fundamental plan e derived for elliptical galaxies 
ins et al.ll2007al [bl). A statistically equivalent for- 
mulation of the BHFP is given by the binding en- 
ergy correlation, which relates the total binding en- 
ergy of the bulge {E ^ M^ a'^) to the mass of t he th e 
BH. Recent observat i ons b y lAUer fc Richstond (|2007D : 
iBarwav &: Kembhavil ()2007[ ) have shown that the bulge 
gravitational binding energy as traced by the stellar light 



profile correlates strongly with the BH mass. 

Figure [T2l shows the binding energy correlation for our 
3:1 and 1:1 disk merger sample (left panel) and for the 
mixed E-Sp and dry E-E mergers (right panel), respec- 
tively, toge ther with the observed relation from the com- 
pilation of lHopkins et al.| ()2007af ). We perform a least- 
squares fit to the data using the equation below 

log(MBH/M0) - e + / log(Af,aV(Moa2)), (6) 

where e and / are the normalization coefficient and slope, 
respectively, and Mq = IQ^^Mq and uo = 200 km/s are 
the normalization values. The best values of e and /, 
together with the dispersion for the various subsamples 
are given in Tab. [TOl We do not find any statistically 
significant difference between the 3:1 and 1:1 merger sam- 
ples. For both samples the derived normalizations and 
slopes are very similar and in good agreement with the 
observed relation. Varying the initial gas mass fraction 
produces more pronounced differences, with a systemati- 
cally higher normalization coefficient for increasing initial 
gas mass fractions. However, the slope remains constant 
within the error bars with changing initial gas mass frac- 
tions. The mixed E-Sp and dry E-E merger samples gen- 
erally maintain the slope of the progenitor sample, with 
a slight increase in the normalization. Finally, combining 
all the merger samples together produces a final relation 
in reasonably good agreement with the observed relation 
and with a similar dispersion to the observed sample. 
The tightness of the derived correlation indicates that 
the final BH mass in all the merger samples is primarily 
set by the depth of the gravitational potential well. 

We conclude that the bulge binding energy relation 
provides a robust predictor for the BH mass over a wide 
range of progenitor masses and initial gas fractions for 
both unequal- and equal-mass disk mergers, mixed E-Sp 
mergers and E-E mergers. The dispersion in the bulge 
binding energy correlation is similar in size to the dis- 
persion in the A/bh — cr relation and both these methods 
should be preferred over the Mbh — M ^, relation that pro- 
duces significantly larger scatter (see lAller fc Richstond 
l2007l for a discussion). 



Equal- and unequal-mass mergers with black holes 
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5. SUMMARY AND DISCUSSION 

In this paper we have studied galaxy mergers with BH 
feedback starting with equal- and unequal-mass disk pro- 
genitors, mixed disk-elliptical progenitors and elliptical- 
elliptical progenitors. Extending on detailed previous 
studies of equal-mass disk mergers (e.g. [Robertson et al.l 
l200 6a^ our merger sample includes 18 unequal-mass disk 
mergers, 16 mixed E-Sp mergers and 16 dry E-E mergers 
in addition to 18 equal-mass mergers resulting in a total 
sample of 68 mergers. The simulations include radiative 
cooling, star formation and black holes wi th their asso- 
ciated f eedback processes as described in Spring el et al.l 
(12005b). 

Wc study the termination of star formation in merger 
simulations including BH feedback. We find that the 
termination of star formation by BH feedback is signifi- 
cantly less important for unequal- mass disk mergers com - 
pared to equal-mass disk mergers (IS pringel et al.l[2005ar ). 
The timescale for star formation termination systemat- 
ically increases with increasing progenitor mass ratios. 
Similarly, a systematic increase is seen in the half-mass 
growth timescales of the BHs, with this timescale varying 
from ~ 0.1 Gyr for equal-mass mergers to 1 Gyr for 6:1 
mergers. This systematic trend can be used as input in 
modeling BH accretion more re alistically in semi-an alytic 
galaxy formation models (e.g. ICroton et al.ll2006l ). For 
mass-ratios of 3:1 and higher mergers with BH feed- 
back are unable to completely quench the star forma- 
tion, with the merger remnants showing star formation 
rates roughly on the pre-merger level even 1 Gyr after 
completion of the merger. 

In mixed E-Sp mergers and dry E-E mergers the star 
formation and BH accretion rates are in general low due 
to the low cold gas mass fraction. The star formation 
is efficiently terminated in mixed E-Sp mergers due to 
the presence of the fully grown super-massive BH in the 
E-progenitor, with lower mass ratio mergers (2:1) again 
resulting in shorter termination timescales compared to 
higher mass ratio mixed (4:1) mergers. Finally, dry E-E 
mergers with two super-massive BHs are extremely effi- 
cient in terminating any residual star formation in the 
progenitors. Remergers of 1:1 merger remnants show 
typically no star formation, whereas remergers of lower 
mass ratio merger remnants can show star formation on 
the level of a few solar masses per year that is rapidly 
terminated in the merger. This implies that already a 
second generation of dry merging is expected to show 
very weak, if any, signatures of residual star formation. 

We show that in addition to 1:1 disk mergers 
(|Di Matteo et al.ll2005[ ) also 3:1 unequal-mass disk merg- 
ers, mixed E-Sp mergers and dry E-E mergers reproduce 
the observed Mbh — o" relation, with the dominant source 
of scatter arising from variations in the initial gas mass 
fractions of the galaxies. The 3:1 unequal-mass mergers 
require BH repositioning, in particular at low numerical 
resolution, which ensures that the final BH quickly reset- 
tles in the central gas disk after the merging of the BHs. 
In addition, the normalization of the observed Mbh — Af* 
relation is also well fitted in our simulation sample, but 
with the resulting slope being slightly steeper compared 
to the observed slope. This discrepancy might be re- 
solved by employing a more aggressive supernova-driven 
wind feedback, as in recent self-consistent cosmological 



simulations (jSiiacki et al.ll2007t[DrMatteo et an2008f) or 

by the fact the progenitors of lower mass galaxy mergers 
have larger initial gas mass fractions compared to high 
mass galaxy mergers. Finally, our complete merger sam- 
ple fits well the observed bulge binding energy correla- 
tion with black hole mass, Mbh — M^cr^, which is statis- 
tically equivalent to t he Black Hole Fundamental plane 
(|Hopkins et al.ll2007airb[ ). The tightness of this correla- 
tion indicates that the final BH mass in all the merger 
samples is ultimately set by the depth of the gravitational 
potential well. 

The remarkable robustness of the simulated Mbh — cr, 
Mbh — M, and Mbh — M^cr^ relations for an extremely 
wide range of i nitial g alaxy properties implies that the 
ISpringel et al.l (|2005bl T effective BH feedback model self- 
regulates the mass growth of the BHs in accordance 
with the observed relations for a very wide range of 
merger progenitors. However, in addition it is impor- 
tant to study the detailed properties of merger rem- 
nants including BH feedback. It has been shown that 
the stellar orbits, e.g. the relative population of box 
and tube orbits, in merger remnant s are affected by the 
amount of gas driven to the cen ter (jBarnes fc Hernquis^ 
119961: iJesseit et all l2005l I2007D . Here, even relatively 
small amounts of gas (~ 5%) can influence higher or- 
der photometric a nd kinematical pro perties of the rem- 
nants as a whole (|Naab et al.ll2006a[ ). As efficient feed- 
back from black holes reduces the central gas content 
in first generation disk merger remnants, future mergers 
will be less dissipative and favor the forma tion of triax- 
ial st e llar system s dominated by box orbits (jJesseit et al.l 
20051: iNaaij et a l. 2006a). In addition, we studied in 



Burkert et al.l ( 2007.) the intrinsic anisotropy S and flat- 



tening e of simulated merger remnants with and with- 
out BH feedback. Here, an equally g ood fit to el- 
lipticals observed by SAURON ( EmscUe nTet all l2004l : 
ICappellari et al.ll2006( l was found in both models that in- 
cluded BH feedback and models without BHs. However, 
binary disk mergers and their descendants were shown 
to preferentially form anisotropic systems with relatively 
high ellipticities, not consistent with observed isotropic 
and round ellipticals. This is kinematical evidence that 
some ellipticals cannot have formed by mergers of disks 
in the first place. In addition alternative formation sce- 
narios ar e motiv ated based on stellar population consid- 
erations ()Naab fc Ostrikcr 2007) as well as direct cosmo- 
logical simulations (jNaab et al.ll2007[ ). 

AGNs are emerging as an important ingredient in 
our current galaxy formation models. Research on un- 
derstanding the observed Mbh — tr, Mbh — M* and 
Mbh — M*(T^ relations using numerical techniques should 
therefore be a useful step in our quest to formulate a more 
complete theory of galaxy formation. 
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APPENDIX 

A. BLACK HOLE FEEDBACK MODEL 

The growth of the supermassive black holes is modeled following the lSpringel et aP (|2005bl ) model. In this effective 
subresolution model the unresolve d accretion onto the BH is related to the resolved gas dis tribution using a Bondi- 
Hoyle-Lyttleton parameterization (|Hovle fc Lvttletonlll939t iBondi fc Hovld[l944l : iBondil [19521 ) . In this description, the 
accretion rate onto the BH is given by 

^""-^ciTv^- ^^^^ 

Here p and Cg are the density and sound speed of the surrounding gas, respectively, v is the velocity of the BH relative 
to the surrounding gas and a is a dimensionless parameter setting the efficiency of the accretion. In addition, it is 
assumed that the maximum accretion is limited to the Eddington rate 

M.,, . ^JI^M^^ (A2) 

where is the proton mass, ctx is the Thomson cr oss-section and Cr is the radi ative efficiency, that we assume to be 
0.1 based on the mean value for radiatively efficient IShakura fc SvunvaevI (119731 ) accretion onto a Schwarzschild BH. 

Thus the accretion rate is given by Mbh = min(MB, Medd)- Finally some fraction ef of the radiated luminosity 
couples thermally and isotropically to the surrounding gas in the form of black hole feedback energy resulting in 

-Efccd = e/^r = e/EiMBHC^- (A3) 

Again, following iSpringel et al.l (|2005bD . we fix ef — 0.05 so that in total 0.5 per cent of the accreted rest mass energy 
is available as feedback. This choice fixes the normalization of the Mbh ~ cr relation and repro duces the o bserved 
relation in equal-mass binary galaxy merger simulations with BH feedback as was shown bv Di Matteo et al.l p005i) . 
The feedback energy is distributed as thermal energy to the surrounding gas particles weighted with the SPH kernel. 
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Fig. B1. — Comparison of the total star formation rate (top), the total bl ack hole accretion rate (middle) and the evolution of the total 
black hole mass (bottom) as a function of time for a 1:1 merger run using the lSpringel et ahl Il2005bl) BH feedback S05 model and the model 
presented in this paper (JOS model). The filled circles indicate the time of merging of the BHs. 



Numerically the BHs are represented as collisionless 'sink' particles, which only feel gravitational forces. The properties, 
such as density, temperature and the bulk velocity of the local gas around the BH are estimated in a similar fashion 
to normal SPH particles. Based on these quantities and the equations above, the BH accretion rate is then estimated. 

B. COMPARISON OF THE BH FEEDBACK MODELS 

We perform a direct compa rison of our BH feedback model to the lSpringel et all ()2005bD model by running one of the 
initial conditions used in the ISpringel et al.l (|2005b r) study (kindly provided by Volker Springel) using our simulation 
code. We use identical numerical parameters governing the stellar feedback, BH feedback and integration accuracy. In 
Fig. IB II we plot the resu l ting st ar formation rate, BH accretion rate and total BH mass for the simulation performed 
using the ISpringel et al.l (|2005bf ) (SOS model, data provided by Volker Springel) and compare it to the corresponding 
output of the model presented in this paper (JOS model). 

The agreement between the two codes is very good, especially given the fact that both the star formation and 
BH accretion processes are treated stochastically. The star formation histories agree very well up to the time of the 
merging of the BHs. After the merging of the BHs we see small differences in the decay rate of star formation in the 
two models. This slight discrepancy is due to the differences in the corresponding BH accretion rates resulting in the 
final BH being marginally more massive in the S05 model. The late decay of star formation is relatively sensitive to 
the final BH mass, with more massive BHs resulting in more efficient termination of star formation and hence small 
differences in the late star formation history are not unexpected. 

We thus conclude given the stochastic nature of the modeled processes that the resulting agreement on the 10% level 
between the two codes is very good. Thus we e xpect our BH feedback prescription to produce nearly identical results 
compared to the original [Springel et al.l (|2005bf ) model. 
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Fig. CI. — 3:1 unequal-mass merger simulations run at the standard numerical resolution with the standard BH merger prescription 
(black lines) and the repositioning method (red lines). The top left and right panels shows the relative distances and velocities between the 
13Hs as a function of time. The bottom panels show the evolution of the BH masses and accretion rates as a function of time. The solid 
line depicts the primary BH, the dashed line the secondary BH and the thick solid line the merged BH, with the time of merging indicated 
with solid circles. 

C. COMPARISON OF THE MERGER PRESCRIPTIONS 

We plot in Fig. [CT] the black hole mass growth and mass accretion rate together with the distance and relative 
velocity between the black holes for one unequal-mass 3:1 merger at our standard resolution using both the standard 
and repositioning merger prescriptions. The galaxies go through their first and second passages at ti = 0.6 Gyr and 
t2 = 1.4 Gyr respectively. The galaxies then approach each other for a third time at ts = 1.8 Gyr at which time the 
BHs in the simulation with BH repositioning merge, unlike the standard case in which the BHs do not merge. The 
mass growth of both the primary (solid lines) and the secondary (dashed lines) BHs are similar for both methods until 
the merger, after which time essentially only the repositioned BH grows in mass. A similar result can be seen in the 
BH accretion rates, which are similar until ~ 1.8 Gyr, after which the repositioned BH shows almost an order of 
magnitude higher accretion rates compared to the standard BH merger case. One can also see that the secondary 
BH mass accretion rate dramatically drops at this time, hence terminating the growth of the secondary BH. This is 
caused by the fact that the secondary BH is perturbed away from its surrounding gas disk during the final stages of 
the merger. 

In the standard BH merger prescription the merging of the BHs is typically delayed with respect to the BHs in 
the repositioning runs. During the BH merger the momentum is conserved, which can give the BH a relatively large 
'kick', especially in the 3:1 mergers where the mass ratio of the BHs is relatively large. The repositioning ensures that 
after the merger the BH quickly resettles in the center of the central gas disk that is feeding the BH. In the standard 
prescription the merged BH remains longer displaced from the central gas disk, thus lowering the mass accretion rate, 
which scales linearly with density (Mb oc p in Eq. lAip . The repositioning also lowers the relative velocity of the BH 
with respect to the surrounding gas. For the standard merging case the relative velocity can be relatively large with 
the BH plunging through the central gaseous disk. This also affects the BH accretion rate as it scales as the inverse 
third power of the relative BH-gas velocity for a fixed sound speed of the medium (A/b v^'^ in Eq. lAip . Finally the 
BH accretion rate scales with the square of the BH mass (A/b oc in Eq. lAip resulting in higher accretion rates 

for more massive BHs, thus an earlier merging of the BHs produces typically larger final BH masses. 

In Fig. IC2I we plot the corresponding details for a 1:1 merger at our standard numerical resolution using the 
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Fig. C2. — 1:1 equal-mass merger simulations run at the standard numerical resolution with the standard BH merger prescription (black 
lines) and the repositioning method (red lines). The top left and right panels shows the relative distances and velocities between the BHs 
as a function of time. The bottom panels show the evolution of the BH masses and accretion rates as a function of time. The solid line 
depicts the primary BH, the dashed line the secondary BH and the thick solid line the merged BH, with the time of merging indicated 
with solid circles. 

two merging prescriptions. The differences between tlie merger prescriptions compared to the 3:1 merger are less 
pronounced with the BHs merging shghtly earher in the repositioned simulation. However, the growth of the black 
hole masses show differences with the BHs in the repositioned run showing roughly equal mass growth prior to the BH 
merging compared to the standard case, where the more massive BH is an order of magnitude more massive than the 
secondary BH at the time of the BH merger. In the standard merging case the secondary BH is perturbed away from 
its gaseous central disk prior to the final merger, thus lowering the accretion rate and growth of the BH in comparison 
to the primary BH. After the merger the repositioned BH exhibits strong growth in a gas-rich environment until the 
thermal feedback energy input is large enough to clear away the remaining gas, thus terminating the BH growth. 
This is not mirrored in the standard merging case, where the primary BH grows strongly prior to merging and only 
modestly after the merger due to the fact that the merged BH is unable to resettle in the remaining central gas disk. 
However, the end result of both merging prescriptions is a final BH mass and stellar velocity dispersion that agree 
with each other on the 10% level. 

Taken together these results indicate that repositioning the BHs produce converged and stable results for the mass 
growth and accretion rates of the BHs. Repositioning is essential in representing the evolution of BHs in 3:1 mergers 
and also produces converged results for 1:1 mergers even at relatively low numerical resolution. 



